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ABSTRACT 

I Aims. We calibrate the number density, completeness, reliability and the lower mass limit of galaxy-cluster detections through their thermal 
SZ signal, and compare them to X-ray cluster detections. 

Methods. We simulate maps of the thermal SZ effect and the X-ray emission from light cones constructed in a large, hydrodynamical, 
cosmological simulation volume, including realistic noise contributions. The maps are convolved with linear, optimised, single- and multi-band 

■ filters to identify local peaks and their signal-to-noise ratios. The resulting peak catalogues are then compared to the halo population in the 
' simulation volume to identify true and spurious detections. 

, Results. Multi-band filtering improves the statistics of SZ cluster detections considerably compared to single-band filtering. Observations 
with the characteristics of ACT detect clusters with masses M > 6 - 9 x 10'^ Mq/H, quite independent of redshift, reach 50% completeness 
at ~ 10'"* Mq/A and 100% completeness at ~ 2 x lO'"* M^/h. Samples are contaminated by only a few per cent spurious detections. This is 

■ broadly comparable to X-ray cluster detections with XMM-Newton with 100 ks exposure time in the soft band, except that the mass limit for 
' X-ray detections increases much more steeply with redshift than for SZ detections. A comparison of true and filtered signals in the SZ and 
, X-ray maps confirms that the filters introduce at most a negligible bias. 



Key words. Cosmology: theory, cosmic microwave background - Galaxies: clusters: general - X-rays: galaxies - Methods: N-body simulations 



1. Introduction 



Galaxy clust ers efficiently trace the large-scale structure of the Lope z-Caniego et"al] 2006 



20061: iHerranz et alJ l2002allbh and investigated the statisti 
cal properties of SZ detect i ons in simu lated m aps (see, e.g., 



Universe (Eisenstein et alj | 2005; lijitsi 2006), and their mass ,2006c ISchafer & Bartelmannll2006l) . These studies were based 



function strongly depends on the cosmological parameters 
dWel er & Battvel 120031: ISefusatti et al.ll2007l: Ipang & HaimanI 



Mehn et all 120061: IVale & White 



20071). They can be easily studied at different frequencies, 



ranging from microw aves to X-rays, and via weak and strong 
gravitational lensing (Clowe et al. '2004'; 'Hen nawi & Spergel 
[2005; Tang & Fan 2005; Co merford et al. 200 6). The Sunyaev- 
Zel'dovich (hereafter SZ) effect and the X-ray emission are 
complementary, and combining them permits individually im- 
possible constraints. Moreover, the SZ effect allows cluster de- 
tections to larger distances than X-ray emission. In turn, grav- 
itational lensing probes the total mass rather than the gas, and 
is most sensitive at intermediate distances between the sources 
and the observer. 

Expecting wide-field data in the near future, many authors 
developed different t echniques to detect clusters through their 
thermal SZ signal dDiego et alj 120021 : iLopez-Caniego et al 



on analytic or semi-analytic cluster models applied to cos- 
mological simulations, but gas physics was not included, 
which is known to affect the appearance of clusters (see, e.g. 



Kazantzidis et al. 2004; Puchwein et al. 2005; Bonaldi et al 



20071). Physical and statistical properti es of X-ray galaxy 



clusters have been extensively studied ([Bprgani et alJ 12004 



RoncareUi et al.1 120061: iMazzotta et al.ll20()4l) . showing that X- 



ray cluster properties are very sensitive to gas physics (see, e.g. 



Carlstrom et"ani2002l) . 



In this paper, we investigate the statistical properties of 
cluster detections in simulated maps via single- and multi-band 
matched filters, applied to X-ray and thermal SZ maps. We 
study completeness and contamination of cluster catalogues 
and the lowest detectable cluster mass. We also cross-correlate 
X-ray and SZ cluster detections. The maps are simulated by 
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projecti ng the outputs of a c osmological hydrodynamical sim- 
ulation ( Borgani et alJ2004 ) of the concordance ACDM model 
along deep light cones. The simulation contains gravitation and 
gas dynamics and includes also several physical processes that 
affect the baryonic component, such as a star formation, su- 
pernova feedback, radiative cooling and a simple reionisation 
scenario. 

Our final goal is to optimise cluster detections in terms of 
reliability and completeness with linear filtering techniques ap- 
plied to multi-wavelength data. Such cluster catalogues provide 
information on the underlying cosmology, the statistical prop- 
erties of the cluster population and its redshift evolution, and 
thus on the non-linear structure growth in the universe. 

We summarise the SZ effect, the X-ray emission and our 
models for them in Sect. |2] The cosmological simulation is 
briefly described in Sect.[3]together with our procedure for cre- 
ating SZ and X-ray maps including noise. The single and multi- 
band filters are described in Sect. |4] In Sect. |5] we analyse the 
properties of the synthetic catalogues obtained by applying the 
single- and multi-band filters, while the correlation between SZ 
and X-ray detections is discussed in Sect. |6] Our main results 
are summarised in Sect.]?] 

2. The observables 

2.1. The SZ effect 

The SZ effect is due to the inverse Compton scattering of CMB 
photons off' the intracluster electrons. The thermal SZ (tSZ) ef- 
fect is caused by the thermal motion of the electrons, the ki- 
netic SZ (kSZ) effect by the bulk motion of the clusters. The 
kSZ effect is typically an order of magnitude weaker than the 
tSZ effect. The tSZ effect is quantified by the dimension-less 
Compton y-parameter 

ksCTT 



y = 



J" dlrie 



(1) 



integrating over the product of the electron number density tig 
and their temperature T,,, where ks is the Boltzmann constant, 
cTj the Thompson cross-section, m,, the electron rest mass and c 
the speed of light. The tSZ effect changes the CMB brightness 
temperature by 

AT 

= ygvix) , (2) 



where gv(x) is the tSZ spectrum 



gv(x) = 



-4 



(3) 



tanh(A:/2) 

as a function of the dimension-less CMB photon energy x = 
hv/iksTcMB)- In the Rayleigh-Jeans limit (x <K 1) gvix) ^ -2, 
and the tSZ effect vanishes at v ^ 217 GHz. 

The kSZ effect causes brightness-temperature fluctuations 

AT 



I alrieVr 

c J 



(4) 



where v, is the radial component of the cluster velocity, defined 
positive for receding clusters. The kSZ effect does not depend 
on frequency and is thus best estimated from observations near 
V ^ 217 GHz, where the tSZ effect vanishes. 



2.2. The X-ray emission 

X-ray emission of galaxy clusters due to thermal free-free 
emission (bremsstrahlung) has the bolometric emissivity 

1.2 X IQ-^'^ T^'\lg(T) 



(5) 



in cgs units, where the gas temperature T is given in keV, np is 
the proton density in cgs units, and g{T) is the Gaunt factor The 



thermal free-free emission spectrum is /(v) = e 



-hv/ksT 



giv, T), 



where y - Z/Zq is a pa rameter depending on the metalhc- 
ity (Z a; 0.3 Zq) (see e.g. iRosati et al.ll2004i Hashimoto et al. 
2004). 

The X-ray luminosity in the energy interval [£1 ,£2] is given 
by Lband = ' ^tonf ^^cre L^°' = egV, and V is the emitting 
volume. The band-function is defined as 



plE,.E: 
band 



■■«2.o(l+Z) 



/^■«2.0U + 

-^.fi.o(i+z; 



C fix) dx 



(6) 



where x is defined as above and C is a normalisation constant 
defined by C ^ f(x)dx - 1. The subscript in the integra- 
tion boundaries means that the band limits are taken in the ob- 
server's rest-frame, thus the (l + z) term in Eq. (|6]l represents 
the K-correction. We shall consider three different bands: soft 
[(£i,£2) = (0.5,2) keV], hard [(£i,£2) = (2,4) keV], hardest 
[(£i,£2) = (4,10)keV]. 

3. The simulations 

3.1. The cosmological simulation 

To create realistic mock catalogues of SZ and X-ray halos, 
we u se the outputs of a h ydrodynamical cosmological simula- 
tion ( Borgani et al.l2004 ) carried out using the GADGET-2 code 
(ISpringeL2005i) . Here we briefly report the main characteristics 
of the simulation. The SZ effect and X-ray properties of objects 
identified in this simulation are described in other studies (e.g . 



Diaferio et al.ll2005t iRoncar elh et al 
Murante et al. 2004; Rasia et al.ll2005 



2OO7I: 'Ettori etal 



iRoncarelli et al.l 



2004 



20061) . 



where further detail can be found. 

The simulation assumes a concordance ACDM model with 
matter density parameter Q,„ = 0.3, cosmological constant 
Qa = 0.7, and Hubble parameter h = (//q/IOO km/s/Mpc) - 
0.7. The power spectrum is normalised to cr^ = 0.8. These 
cosmological parameters are consistent with recent observa- 
tional estimates derived from CMB and w eak lensing data (e.g. 
Spergel et al.ll2007tlHoekstra et al.ll2006h . 

The computational box with a comoving side length of 192 
Mpc/h is filled with 480^ particles each of dark matter and gas. 
The mass resolution is thus 4.6 x IO'^Mq/Zi and 6.9 x IQ^M^Ih 
for dark matter and gas, respectively. The Plummer-equivalent 
gravitational softening is set to epi - 7.5 kpc/h comoving 
for < z < 2, and switches to physical units at higher red- 
shifts. The gas component is treated including several physi- 
cal processes: a hybrid m ulti-phase model for star form ation in 
the interstellar medium ( Springel & HernquistI 2003afl . radia- 
tive cooling within an optically thin gas consisting of 76% of 
hydrogen and 24% of helium by mass, supernova feedback to 
model galactic outflows, and heating by a time-dependent, pho- 
toionising uniform UV background given by quasars reionising 
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the Universe at z 6 (iHaardt & Madaulll996l) . The output of 
the simulation consists of one hundred snapshots, logarithmi- 
cally equidistant in redshift between z - 9 and z = 0. We use 
these snapshots to create our mock light-cones. 

3.2. Construction of the halo catalogues 

To identify the halos in each snapshot, we run a friends-of- 
friends algorithm with a linking length of 0.15 times the mean 
particle separation to identify particles belonging to the same 
group. Next, we identify the halo centre as the position of the 
particle at the local minimum of the gravitational potential. The 
final halo catalogue contains the positions, the virial masses, 
radii and redsh ifts of all halos in each snapshot. More details 
can be found in Pace et al. (2007). 
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Fig. 1. Average number of halos per mass bin per square de- 
gree. The solid and dashed lines show the numbers for light- 
cones extending to z = 1 and z = 2, respectively. The number of 
high-mass halos is approximately the same in both case, while 
the number of low-mass halos is higher for the deeper cone. 
The values are averaged over eleven different realisations. 

Figure [T] shows the mass function of the halos contained in 
light cones for limiting redshifts of z = 1 and z = 2, averaged 
over eleven different realisations and normalised to one square 
degree. Light-cones reaching z = 2 contain more low-mass ha- 
los than light-cones up to z = 1 because of their larger volume, 
while the number of massive halos is almost unchanged since 
most of them are at low redshifts. 



3.3. Construction of the simulated SZ and X-rays 
light-cones 

To construct a realistic three-dimensional distribution of mat- 
ter, we construct several light cones by stacking the snapshots 
of our cosmological simulation at different redshifts, observing 
two aspects: First, we shift and rotate each snapshot to make 
it completely independent from the others. Second, we cut the 
boxes appropriately to avoid including matter from more than 
one snapshot, since contiguous snapshots partially overlap in 
redshift. 



To create the maps, we proceed in two steps. We first 
project gas particles on a two-dimensional grid, and then stack 
the resulting planes summing the contribution from each indi- 
vidual plane to obtain the final map. Projecting the gas parti- 
cles on a regular grid, we take the smoothing length of each 
particle into account and use the SPH (Smoothed Particle 
Hydrodynamics, see e.g. Monaahan 1992) inter polating kerne l 
from the original simulation dMonaghan & Lattanzid Il985h . 
Particles are distributed over grid cells with a fractional weight 
calculated from the overlap of the spline kernel with the grid 
cell. 

Since each particle has its individual smoothing length, the 
number of pixels involved in each projection varies, which pro- 
du ces appropriately smoo th maps. This algorithm was applied 
bv lPuchwein et al. I (l2005h to stu dy the impac t of gas physics on 
strong lensing by clusters and bv lda Silva et a l. ( 2000) to create 
simulated SZ maps. 

We repeated this procedure for each snapshot of the simu- 
lation, and then summed the contributions of all planes to ob- 
tain the final map. We created eleven different maps including 
matter up to redshift z = 1 and z = 2 with a resolution of 
2048 X 2048 pixels. The opening angle of the light-cone is de- 
termined by the last plane of the pile, because its dimension is 
fixed in comoving units: it corresponds to 4.9 and 3.1 degrees, 
for light-cones up to z = 1 and z = 2, respectively. 

We built our light-cones in the same way as for 
gravitational-lensing maps, i.e. projecting the whole box on 
a regular grid and then selecting the p ortion of the plan e en- 
closed by the light cone (see Fig. 1 in IPace et alJ (120071) ). In 
other words, our planes all have the same number of pixels 
with constant comoving size, chosen acc ording to the simu- 
lation resolution. Oth er authors (see e.g. Ida Silva et al. I I2OOOI: 
Roncarelli et al ] l2006 l) prefer to construct the light-cones keep- 



ing the angular resolution of the pixels constant, such that each 
plane contributes the same amount of pixels. We avoid this be- 
cause then the pixel resolution would exceed the physical res- 
olution of the cosmological simulation, especially at low red- 
shift. 

3.3.1 . Simulating tine SZ maps 

In order to create a map for the SZ effect we replace the line- 
of-sight integrals in Eqs. ([T) and (HJi by sums over gas particles. 
The contributions to the tSZ and kSZ effects of the ;-th particle 
are thus 



1 kB(TT 

iJ. rrieC^ 

pix 

1 CTj- 



nejTi 



pix 



(7) 



(8) 



respectively, where Lpix is the physical size of a pixel at the 
distance of the corresponding plane. We approximate the radial 
velocity vv,/ with the velocity component along the z-direction 
and relate the temperature of the gas particles T (in Kelvin) to 
their internal energy per unit of mass {V in km^/s^) through 



r = 10" X —mpii U 

iKB 



(9) 
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as for a monatomic ideal gas, where irip is the proton mass and 
/i is the mean molecular weight 



1 + 4yHe 



(10) 



I +yHe+ rie 

We adopt jH^. ~ 0.08. 

3.3.2. Simulating tine X-ray maps 

We create the light-cone maps for the X-ray emission in the 
three energy bands mentioned above: soft (0.5 - 2 keV), hard 
(2-4 keV) and hardest (4 - 10 keV). We model the contri- 
bution of each particle to the X-ray signal in the soft a nd hard 
bands adopting the M eKaL model (|Mewe et al. 1995 ) as im- 
plemented in XSPEC ( )Arnaudl[l996r It returns the emission 
spectrum of the hot diffuse plasma and is particularly suited for 
the soft band where the influence of the metal line emission is 
important. Using this model, the X-ray luminosity is 



J^x,i = {nipH) ^miPiXeA(Ti,Zi, E[, E'2) 



(11) 



where Xe = rielnn is the ratio between the number density of 
free electrons and hydrogen nuclei. The cooling function A de- 
pends on the particle temperature, on the metallicity and on the 
energy band [£'i,£'2]- 

For the hardest b and, we can use the model described in 
Borgani et al. (Il999h because the influence of metal lines is 



negligible there, and a simple power-law parametrisation suf- 
fices. It is then possible to parametrise X-ray luminosity of the 
;-th particle in the hardest band by 



I hardest 



1.7 X 10^ 



erg/s 



(12) 



where m,, p, and T, are its mass, its density, and its temperature 
in keV. 

For each gas particle, the X-ray flux is 



7 I I 2 

" TT7^ erg/s/cm 



(13) 



where di(z) is the luminosity distance of the particle from the 
observer. As for the SZ maps, we project the X-ray flux on a 
regular grid using the SPH kernel. 

3.4. Simulating observations 

We simulate observations with specific instruments, namely the 
Atacama Cosmology Telescope (ACT) for the SZ efifect be- 
cause of its high angular resolution and sensitivity, and the 
space telescopes XMM-Newton and Chandra for the X-ray 
emission, assuming exposure times of 30 and 100 ks. 




100 150 
V [GHz] 



300 



Fig. 2. The spectrum of the tSZ efifect, with the ACT frequency 
bands overplotted (crosses). 



filter because it contains precious information on the noise due 
to the primary CMB anisotropics. 

The noise in tSZ observations has three main contributions: 
the CMB radiation, the instrumental noise, and the kSZ effect. 
The noise due to the primary CMB anisotropics and to the in- 
strument can be modelled as two independent Gaussian random 
fields. We us ed the CMB power spectrum as computed with 
CMBEASY (lDoranll2005b . and the instrumenta noise power 
spectrum 



- w exp 



l{l+l)FWHM'^ 



81n2 



(14) 



where FWHM is the full width at half maximum of the instru- 
mental beam in arcmins, I represents the multipole order, and 
= (AT/T FWHMf. The quantity AT/T gi ves die sensi - 



tivity of the instrument on the scale of the beam ( Knoxl[l995 ). 
Finally, the maps are convolved with the instrumental beam. 
Table [1] summaries some parameters used to mimic ACT ob- 
servations. 



Central Band Frequency 


FWHM 


AT/beam 


(GHz) 


(arcmin) 




145 


1.7 


2 


225 


1.1 


3.3 


265 


0.93 


4.7 



Table 1. Parameters used to produce mock observations with 
ACT. In the first column we show the central bands frequen- 
cies at which ACT operates; second and third columns list the 
corresponding FWHM and sensitivity, respectively. 



3.4.1 . Noise inclusion into the SZ map 

Figure |2] shows the tSZ spectrum gvix). The crosses indicate 
the three frequency bands of ACT. The v - 225 GHz channel 
is close to the frequency where gv{x) vanishes and the dominant 
signal will be due to the kSZ effect. Nonetheless, we shall see 
that this channel is very important for our multi-band frequency 



In Fig. |3] we show an example of the SZ simulations. 
The right panel displays the simulated observation with the 
V - 145 GHz channel of ACT, including the kSZ effect, the 
instrumental noise and the primary CMB fluctuations. The left 
and central panel show instead a simulated map for the tSZ and 
kSZ effect, respectively. 
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Fig. 3. Examples of SZ maps including cosmic structures to z = 1 . The left panel shows the simulated tSZ map, the central panel 
the kSZ map, and the right panel shows a simulated observation of the same field in the 145 GHz channel of ACT, including the 
kSZ effect, the CMB signal and the instrumental noise. The size of the maps is ^ 4.9 deg^. 



3.4.2. Noise inclusion into tine X-ray maps 

In order to transform our X-ray flux maps into count rates, 
we first multiply them by the energy conversion factor of the 
instrument. For each energy band, this factor is computed as- 
suming a spectral model, consisting of thermal free-free emis- 
sion with kgT = 4 keV, convolved with the response matrices 
of the EPIC-PN detector on-board XMM-Newton or with the 
ACIS-I array on-board Chandra. Assuming an exposure time 
(either 30 or 100 ks), we transform our count-rate maps into 
photon-count maps. We then added the background of an X-ray 
observation, given by the sum of the detector noise, the unre- 
solved X-ray background, and the particle background. Finally, 
we convolved the resulting maps with the point-spread function 
(PSF) of the instrument (5 and 0.5 arcsec for XMM-A^ewfo« 
and Chandra, respectively). In order to reproduce the observa- 
tional noise, we added Poissonian noise according to the local 
photon flux. The noise levels used in the simulated observations 
with XMM-Newton and Chandra are summarised in Tab. |2] 




Fig. 4. Examples of X-ray maps. The left panel shows the X- 
ray flux in the soft band for the same field displayed in Fig. [3l 
the right panel represents a simulated observation with XMM- 
Newton assuming an exposure time of 30 ks. The size of the 
image is the same as in Fig. [3] 



Instrument Band Noise Effective area 

(counts/sec/deg^) (atl.4keV) 
XMM-Newton 0.5 - 2 keV 24.624 
(EPIC-PN) 2 - 4 keV 11. - 1 200 cm^ 

2-10 keV 19.64 

Chandra 0.5 - 2 keV 1% 
(ACIS-I) 2-4keV 2.47 -600 cm- 

4-10 keV 7.4 



Table 2. XMM.-Newton and Chandra noise levels and effective 
detector areas for the EPIC and ACIS-I instruments, respec- 
tively. 



In Fig.|4]we give an example of X-ray simulation in the soft 
band. The left panel shows a synthetic X-ray map, while the 



corresponding simulated observation with XMM-Newton, in- 
cluding the background, instrumental and Poissonian noises is 
shown in the right panel. An exposure time of 30 ks is adopted 
here. The maps are based on the same mass distribution under- 
lying Fig. [3] 



4. Filtering method 

We now describe the matched filter used to detect clusters in the 
simulated maps. For the X-ray observations, the noise is not 
correlated across bands or with the noise in SZ observations, 
thus the multi-band filter described above will be applied only 
to the SZ observations. 
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4.1. Single-band filter 



4.3. Filter dependence on the template 



The single-band filter ^ we use is analogous to the optimal 
filter proposed bv lMaturi et al. ( 2005 ). We suppose that the ob- 
servational data D(ff) are given by 



D{0)^S(0) + N(0), 



(15) 



where S (0) - At(0) is the signal from sources whose spatial 
shape is modeled by r, and N(0) represents the noise. The esti- 
mated amplitude of the linearly filtered signal is 



^est(e) 



= J cf-e'D{0'y¥{\0' -0\) 



(16) 



Imposing the conditions that the optimal filter be unbiased 
{{A-est - ^> = 0) and that its variance (cr^ - ((Aest - A)^)) be 
minimal, its shape follows from a variati onal problem with the 
two p reviously described constraints (see lHaehnelt & Tegmark 

mm. 



1 



(271)2 



PNik) 



cP-k 



PN(k) 



(17) 



where the hats denote Fourier transforms. Equation (fTTl i shows 
that the shape of the filter *F is determined by the shape of the 
signal, T, and by the power spectrum of the noise, P^. 

4.2. A simple model for galaxy clusters 

We model the gas distribution for the filter template by the trun- 
cated King profile (King 1962) 



p(x) 



1 



for X < r, 



(18) 



where is the core radius, r, is the truncation radius defined as 
ten times the virial radius r,,, and x - r/r^. Although based on 
the simple assumption that gas follows a static and isothermal 
distribution, this profil e is justified by simulations and X-ray 
obser vations (see e.g. iBorgani et al.l l2004i iFinoguenov et al. 
200lh . 

The tSZ effect (Eq. [TJ is proportional to the density p and 
normalised to the Compton-F parameter 



(19) 



integrated over the solid angle f2. 

The X-ray luminosity is proportional to the square of the 
electron density and thus to the square of the profile (fTST l. It 
is normalised such that the bolometric luminosity follows the 
empirical relation 



Lx,boi = 2.99 X 10''4/i"2 



\6keV) 



erg/s 



(20) 



6keV/ 

( Kitavama & Sutd 1997 ). where T - T^iiirs) is the X-ray tem- 
perature. The luminosity in a given band is computed using 
Eq. (|6]l. We do not assume any explicit redshift dependence. 



Figure|5]shows the profile of the filter used for the SZ observa- 
tions at 145 GHz (upper panels) and for the X-ray observations 
in the soft band (lower panels). In the left panels, we assume a 
template with mass M = lO''* M^/h and illustrate how the filter 
depends on the cluster redshift. In the right panels, we fix the 
redshift of the template to z = 0.5 and show how the profile de- 
pends on cluster mass, based on the two cases 
(solid line) and M = 10'** M^/h (dashed line). The upper right 
panel demonstrates that the shape of the SZ filter is insensitive 
to both the mass and the redshift of the template, while only the 
normalisation changes. 





Fig. 5. The matched filter The left panels show the filter pro- 
files normalized to unity (for a better understanding) for differ- 
ent redshifts for a fixed halo mass of M = IQ^* Mq/H. The right 
panels show the filter profile for different halo masses at a fixed 
redshift of z = 0.5. SZ and X-ray filter profiles are presented in 
the upper and lower panels, respectively. 



This implies that its sensitivity is independent of the par- 
ticular choice of the template. This characteristic results from 
the shape of the noise power spectrum, combined from instru- 
mental noise and the CMB: since their power spectra are inde- 
pendent of both redshift and mass, the minimum of the noise 
(and therefore the maximum of the filter) is always at the same 
wavenumbers. Due to the particular shape of the noise power 
spectrum (shown in Fig.|6]l, the filter needs to oscillate because 
of the lack of information in a single band of the global charac- 
teristics of the background noise. Unlike the SZ filter, the X-ray 
filter (see the lower panels of Fig.|5j is simply proportional to 
the template because the noise is white. This implies that struc- 
tures in the signal do not change. Thus, varying the template, 
the X-ray filter changes accordingly. 

4.4. The multi-band matched optimal filter 

The multi-band matched filter allows combining information 
from differeri t bands. For its complete derivation, we refer to 
I Schafer et"al] ([2006.) and iMelin et all (l2006i) and references 
therein. Generalising the single-band case, the data obtained 
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Fig. 6. The noise power spectra. The sohd Une shows the CMB 
power spectrum, the dashed hne the instrumental noise, and the 
dotted-dashed hne the inverse of the total noise, given by the 
sum of CMB and instrumental noise. The filter is proportional 
to the inverse of the total power spectrum. 

at a given frequency v are the sum of the signal S y{0) and the 
noise A^v(^), 



Dy{0) ^ S y{0) + Ny{0) . 

The signal is modelled as 
S,(0)^Af,T,{0) , 



(21) 



(22) 



where A is the (band-independent) amplitude, fy is the fre- 
quency dependence of the amplitude and Ty{0) is the spatial 
profile (normalised to unity). The index v runs from 1 to the 
number of available bands, M. We assume that the background 
noise has zero mean in each band and that its statistical proper- 
ties are fully characterised by the correlation functions 

-- {NAk)NyJkr} = (2nf6(k - k')P,, 



C 



(23) 



where Pi,v,v2(k) is the cross-power spectrum. 

To measure the signal amplitude A, we define a linear esti- 
mator for each band. 



Aest,v(e) = j d^6'DA0')'i'y(\0' 
and a total estimate by 

M 

''lest — ^e.st,v(^) ■ 
v=l 



•01) 



(24) 



(25) 



In analogy with the single-band filter, 4'(0) = [4'y(0)] is the 
optimal filter, minimising the estimated variance, cr^ = ((Aest - 
<A))^), and avoiding bias, b = (Agst - A) - 0. 
In deriving the filter, it is convenient to define the filter vector 
^'(0) = [*I'^,(6')] and the signal vector F(6') = [Fy(0)]. We further 



de fine the cross-p ower matrix C*'' = [Cy'^^j. 



Schafer et alj 



conditions is 



show that the filter satisfying these 



(26) 



where C ' is the inverse of the matrix C, and the normalisation 
factor a is given by 

a"' = F^C^'F . (27) 



The estimate has the variance 



(2;r)2 



(28) 



According to Eq. ( |26] ), the matched filter depends on the 
noise, taking advantage of its correlation between all bands. In 
the special case Cv,v, = Cy, <5(vi -V2), i.e. in the case of uncorre- 
lated noise, the matrix C and thus also its inverse are diagonal, 
and the optimal filter for each band had the same shape as the 
single-band filter defined in Section l4~n Just the normalisation 
differes. 

Figure|7]shows the matched-filter profile for the three bands 
of ACT and a template with M = lO'"* Ms/h at z = 0.4. It 
is evidently broader than for a single band (see Fig. |5]l. Note 
also that the multi-band SZ filter does not oscillate because the 
multi-band filter has complete information on how to separate 
the signal from the noise pattern and therefore does not need 
to compensate th e background by os cillations. The matched 
filter described in lSchafer et al.l (120061) shows oscillations (see 
Fig. 6) because a different noise power spectrum was used to 
define the filter, in particular Galactic foregrounds were added. 
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Fig. 7. The profile for the matched tSZ filter for a halo template 
with M - lO^'* M0//1 at z = 0.4. The sohd, dashed and dash- 
dotted lines refer to the 145 GHz, 225 GHz and 265 GHz bands, 
respectively. 



5. Results: Statistics of tlie detections 

We now apply the single- and multi-band filters to our set of 
eleven SZ and X-ray simulated maps and statistically anal- 
yse the samples obtained. We define a detection as a group 
of neighbouring pixels in the signal-to-noise (hereafter S/N) 
maps, whose values exceed a certain threshold. The pixel with 
the largest S/N value defines its position. A detection is true 
if it can be associated to a halo present in the original ACDM 
simulation, while it is spurious if caused by noise. 

In order to reliably associate each detection with a halo 
from t he parent simulatio n, we apply the method described and 
used in lPace et al.l (l2007h . to which we refer for details. Briefly, 
we produce sequences of maps by removing sequentially one 
plane from the stack. If a halo is in fact located on an individ- 
ual plane, its detection must disappear from the map when this 
plane is removed. If a detectable peak keeps appearing in all 
maps of a sequence, it must originate from the noise. 
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5.1. SZ single-band detections 

Due to the oscillations of the single-band filter (see the upper 
panels of Fig. |5]l, detections consist of positive S/N peaks sur- 
rounded by negative and positive ring-like structures. Typically, 
these rings are fragmented by noise and appear as secondary 
S/N peaks suiTounding the most prominent detections. These 
secondary peaks can easily be confused with true detections. 
Some examples can be seen in the left panel of Fig. [8] 
Consequently, we need to include these correlated noise struc- 
tures in our noise model. 




Fig. 8. Maps for the S/N ratio illustrating the procedure to sup- 
press artificial structures in the filtered SZ maps using the cor- 
related noise. The left panel shows the detections with S /N >2 
as obtained with the filter, while the right panel shows the noise 
map including the damped oscillations of the single-band filter. 




Fig. 9. S/N map, as shown in Fig. [8] after removal of the spuri- 
ous ring artifacts. The side length of the map is 2.5 degrees. 

This is possible because we can easily estimate their ex- 
pected shapes and amplitudes, which derive from the convolu- 
tion of the signal with the oscillatory pattern of the filter (see 
Fig. |5]l. The noise map is thus obtained adopting an iterative 
approach. First, we compute the S/N maps by applying the fil- 
ter and estimating the noise. Second, we identify the detection 
having the highest S/N value, compute the correlated noise pat- 
tern as explained above and include it into the noise map. Then, 



we use this new noise estimate to compute an updated S/N map, 
where we identify new detections searching for local S/N max- 
ima. We can iterate this procedure until no more detection is 
found with S/N exceeding a given threshold. The noise map re- 
sulting from this procedure is shown in the right panel of Fig. [8] 
The resulting S/N map obtained using the correlated noise pat- 
tern is shown in Fig.|9] where, thank to the use of the correlated 
noise, we notice that the fragmented rings are efficiently sup- 
pressed. 
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Fig. 10. The total number of detections per square degree (up- 
per panels) and the fraction of spurious detections (bottom pan- 
els) as functions of the S /N ratio. The results are averaged over 
eleven realisations and shown for light cones extending to red- 
shifts z - I (left panels) and z - 2 (right panels). Different line 
styles refer to the different bands considered, as indicated in the 
plots. 

The upper panels of Fig. [TO] show the total number of de- 
tections per square degree in our ACT simulations. Results are 
shown for all three channels and for light cones limited to z = 1 
and z - 2 (left and right panels, respectively). The coiTespond- 
ing fractions of spurious detections are displayed in the bot- 
tom panels. Results are averaged over eleven maps covering 
24.5 deg^ and 9.5 deg^ and renormalized to 1 deg^ for limiting 
redshifts z - I and z = 2, respectively. As expected, we find 
very few detections at v = 225 GHz where the tSZ effect is 
very weak. 

Since the SZ effects depend only mildly on redshift, we find 
more objects in a light-cone extending up to redshift z = 2 than 
for z - 1. In particular, at v = 145 GHz and v - 265 GHz, 
we find at least 30% more detections in the deeper light cone. 
The main differences for the results in these two bands are at 
low signal-to-noise ratios, where more structures are detected 
in the highest-frequency band, owed to the higher sensitivity 
of the instrument. For S/N > 5, the number of detections is 
similar in the two bands. 

Because of the low tSZ signal, the fraction of spurious de- 
tections at V = 225 GHz is very high 80% for S /N ~ 3) and 
virtually independent of the limiting redshift of the light-cones. 
For the maps at v = 145 GHz and v = 265 GHz at z = 1 (z = 2), 
we find that 20% and 40% (30% and 40%) of the detections at 
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the higher completeness in the former band. This is due to the 
fact that the band at 265 GHz has a better FWHM (see Table 
[T]) and the noise is therefore lowe r (see E q.fT4l). 



Fig. 11. Upper panels: the fraction of detected halos as a func- 
tion of the halo mass. Lower panels: the sensitivity of the 
method expressed in terms of the minimum detectable mass as 
a function of redshift. Results are shown for light-cones extend- 
ing to z = 1 and z = 2 in the left and right panels, respectively. 
Different line styles refer to the three different ACT bands, as 
labeled. The curves represent average over eleven realisations. 



S /N ~ 3 are spurious, respectively. For S /N > 4, the fraction 
of spurious detections drops to a few per cent. Including matter 
up to redshift z = 2 into the light cones slightly increases the 
fraction of spurious detections because unresolved structures 
overlap. The number of detections with high S /N does not in- 
crease with the depth of the light cones because the massive 
structures responsible for them are mainly located at low red- 
shifts, but the number of low-5 /N detections is increased by up 
to a factor of five for the least massive detectable halos. 

The upper panels of Fig. [TT] show the fraction of halos de- 
tected in different bands as a function of mass. The left and 
right panels show the results for light cones up to z = 1 and 
z = 2, respectively. The values shown are obtained after av- 
eraging over eleven realisations. The lower panels of Fig. [TT] 
present the sensitivity of the detections, expressed in terms of 
the minimum detectable mass as a function of redshift. To re- 
duce the noise, this was computed combining two subsequent 
planes in the light cone and averaging the mass of the ten least 
massive halos detected in the corresponding redshift interval. 

Since the tSZ effect depends only mildly on the source red- 
shift, the completeness for detections up to z = 1 and z = 2 
is essentially the same. The two channels at 145 GHz and 
265 GHz perform equally well, while the completeness is much 
lower at 225 GHz, as expected. The minimum detected mass is 
M ^ 7 X IO'-'Mq//! and the completeness reaches 100% for 
masses M ^ 3 x IO'^^Mq/Zz. The curve is quite steep, showing 
that the completeness drops rapidly as the mass decreases. For 
halo masses > M ^ 1.5 x IO'^^Mq//;, the completeness has al- 
ready dropped to ~ 50%. The minimum detectable halo mass is 
essentially constant for all z, reflecting the near-independence 
of the tSZ effect on distance. Again, the method is much less 
sensitive at 225 GHz than at the two other frequencies. In addi- 
tion, the filter tends to find halos with lower mass at 265 GHz 
than at 145 GHz (blue and red curves, respectively), explaining 



Similar results were found bv lLope z-Canieg o et al. (l2006l) 



using the same approach to detect clusters in simulated Planck 
data. For low detection l imits (corresponding to low S /N ra- 
tios), |Logez:£anie£oetal. (2006) found a very high fraction 
of spurious detections, decreasing quite rapidly towards higher 
flux limits. However, results can only be compared at a qual- 
itative level because the instrument and the band frequen- 
cies differ greatly. In partic ular, only two channels used by 



Lopez-Canie^o et al. I (1200^ have a frequency range close to 



those of ACT. Moreover, we study S /N ratios rather than fluxes 
and account for the correlated ring patterns due to the oscilla- 
tory filter behaviour in the noise model. 

5.2. Statistics of SZ multi-band detections 

Figure [12] shows an example of the S/N maps obtained by ap- 
plying each component of the multi-band filter to the three 
channels at 145 GHz (upper left panel), 225 GHz (upper right 
panel) and 265 GHz (bottom left panel). The final result is ob- 
tained by combining them using Eq. (125] ) to give the S/N map 
of the multi-band detection in the bottom right panel. 




] [ 





Fig. 12. Example of the application of the multi-band matched 
filter to a S/N map. The starting input is the noisy map in the 
right panel of Fig. [3] Different panels refer to the maps filtered 
at different frequencies: 145 GHz (upper left panel), 225 GHz 
(upper right panel), 265 GHz (lower left panel). The final S/N 
map obtained by summing the previous three maps is shown in 
the lower right panel. 

Compared to the results obtained with the single-band fil- 
ter, the S/N ratios are now enhanced by at least a factor of 
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three. As be fore, the S/N peaks in the maps were identified with 
SExtractor (IBertin & ArnoutsI 119961) which allows easy de- 
blending of nearby detections. The number counts per square 
degree as a function of the S/N threshold are shown in the up- 
per left panel of Fig. [T3] The solid and dashed lines refer to 
light-cones extending up to z = 1 and z-2, respectively. 

The lower number of total detections respect to the single- 
band filter at 265 GHz, has not to be interpreted as a signal 
of a poorer performance of the multi-band filter, but it has to 
be ascribed partly to the fact that the multi-band filter is wider 
and increases the confusion which is not entirely removed by 
SExtractor, but most importantly because of the lower sam- 
ple contamination. In fact, the fraction of spurious detections 
is much lower for the multi-band than for the single-band filter 
We find that the number of real, multi-band detections is higher 
by almost one order of magnitude for S /N ~ 3 and by a fac- 
tor of five for S /N ~ 5 than the single-band detections. These 
results do not depend on the depth of the light cones. 




Fig. 13. Upper left panel: the total number of detections per 
square degree as a function of the S /N ratio, obtained with the 
multi-band filter Upper right panel: the fraction of spurious 
detections as a function of theS /N ratio. Bottom left panel: the 
fraction of halos detected as a function of mass. Bottom-right 
panel: the minimum detectable mass as a function of redshift. 
Solid and dotted curves show results for light cones extending 
to z = 1 and z-2, respectively. 

The top-right panel of Fig. [13] shows the fraction of spuri- 
ous detections in the final catalogue. Again, the great advantage 
of multi-band observations is obvious. While the fraction of 
spurious detections in the 265 GHz channel alone is ~ 40% at 
S /N = 3 and ~ 5% atS /N ~ 5, multi-band filtering reduces the 
fraction of spurious detection to the level of a few percent even 
at very low S/N ratios. Thus, the vast majority of the detections 
in the top-right panel is due to real structures. The completeness 
of the catalogue is shown in the bottom left panel of Fig. [13] We 
find that 20%, 50% and 80% of the halos with masses exceed- 
ing 7 X 10'3 Ms/h, ~ 10^'^ Ms/h and 2 x 10^"^ M^/h contained 
in the light cones are detected by the multi-band filter, respec- 
tively. In contrast, the single-band filter detects only ~ 20% 
of the halos with masses exceeding 10^'* M^/h in the 145 GHz 



and in the 265 GHz channels. Finally, we present in the bottom 
right panel of Fig. [13] the minimum halo mass detectable for 
the multi-band filter as a function of redshift. Again, the curves 
show a very weak dependence on redshift, although the val- 
ues seem to weakly increase with redshift. The mass threshold 
ranges from ~ 6 x 10'^ Mq//i at z = to ~ 9 x 10'^ M^/h at 
z - 1.8, which confirms that the multi-band filter also improves 
the sensitivity. 

We conclude this section by comparin g our result s to an 
independent study of multi-band filtering. Melin et al.l (120061) 
adopted a multi-band matched filter for analysing Monte-Carlo 
simulations. Mimicking SZ observations with the South-Pole- 
Telescope (SPT), they performed a statistical study of the de- 
tections comparable to ours. Since the SPT bands and FHWM 
are similar to those of ACT, a fair comparison between our and 
their results is possible. They found 6 detections per square 
degree atS /N > 5, very close to the 8 detections per square de- 
gree we find. The lowest masses they can detect is three times 
higher than we find for ACT. This is expected because of the 
higher noise level of SPT compared to ACT. The co ntamination 



by spurious detections in the catalogues found by iMeUn et al 
([2006,1 is also similar to ours. 



5.3. Statistics of X-ray detections 

We now turn to the halo detections in the simulated X-ray 
maps. We show results only for light-cones extending up to 
z - 1, because the number of X-ray sources at higher red- 
shifts is very low, and the corresponding plots are very similar 
to those for z = 1 . 

The left panels of Fig. [14] show the total number of de- 
tections per square degree. From top to bottom, results are 
shown for the soft, hard and hardest bands. They represent 
the average over eleven different map realisations and numbers 
are normalized to 1 deg^. The number of detections found for 
XMM-Newton is at least five times higher than with Chandra. 
Although XMM-Newton has a higher background level than 
Chandra, its better performance is due to the different effective 
area of the instrument (see Table O. The number of detections 
decreases quite rapidly with the S/N threshold. We find few 
objects with high signal-to-noise ratios also in the high-energy 
band. 

The fraction of spurious detections is presented in the right 
panels of Fig. [14] We notice that it is always below 20% in all 
bands for detections with S /N > 3, and that it drops to a con- 
stant value a: 5% independent of exposure time and instrument. 
This is due to few spurious peaks having quite high S /N ratio. 
Since also the number of real objects with high S/N ratio is ap- 
proximately constant over a large range of S/N, their ratio is 
almost constant. These few, but prominent spurious detections 
are due to regions in the maps where several low-mass halos 
are approximately aligned along the line-of-sight. 

Figure [15] shows the completeness of the detections. In the 
soft and hard bands, the detection catalogues for XMM-Newton 
are always more complete than for Chandra. In the soft band, 
the completeness decreases much more slowly as a function 
of mass than in the harder bands. The completeness reaches 
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Fig. 14. Total number of X-ray detections per square degree 
(left panels) and fraction of spurious detections (right panels) 
as a function of the S /N ratio. The results are shown for three 
different bands analysed: soft band (upper panel), hard band 
(middle panel), hardest band (bottom panel). The light cones 
extend to z - 1. Different line styles correspond to Chandra 
and XMM-Newton, respectively), different colours to different 
exposure times (black and red curves are for 30 ks and 100 ks, 
respectively). The results are averaged over eleven reahsations. 



100% for masses M > 2 x 10^* Me/h, and it is ~ 50% at 
M ~ 3 X 10'^ Mq/H for XMM-Newton with an exposure time 
of 100 ks, while it is of the order of few percent for Chandra. 
In the hard band, Chandra detects only halos with mass M > 
2- 3 X 10'"* Mq/Zz, and the detection threshold drops by a factor 
four with XMM-Newton. Results for the soft band are similar, 
but the softer photons of lower-mass halos lower the detection 
limit to M a! 10'^ Mq/It. Conversely, only halos with masses 
exceeding a few times lO'"^ Mq/H can be detected in the hardest 
band. There, the completeness for Chandra is almost as high 
or even larger than for XMM-Newton for an integration time 
of 100 ks, because only very massive halos are detected, and 
Chandra has a better resolution than XMM-Newton (see Table 
|2|l. As the X-ray emission depends mostly on the particles of 
the halo core, it might be that X-ray properties for halos of 
mass M x IO'^Mq/Zi are not very well converged, therefore 
reliable results can only be safely inferred to halos of mass M ^ 
5 X 10'^ Mq/Iz. The limit is shown in the left panels of Fig.fTSl 
by the black vertical solid line. 

The right panels of Fig. [15] show the sensitivity for the X- 
ray detections. In the soft band, Chandra performs worse than 
XMM-Newton and thus requires larger masses for detection. 
In the hard and hardest bands, the sensitivity at low redshifts 
is similar for both satellites, even if Chandra performs slightly 
better because of its better resolution. Chandra finds halos only 
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Fig. 15. Fraction of halos detected as a function of halo mass 
(left panels) and sensitivity of the method (expressed in terms 
of minimum detectable mass, right panels) in three different 
X-ray bands: soft band (upper panel), hard (middle panel), and 
hardest (bottom panel). Light cones extend to z = 1 . Line-styles 
and colours refer to different instruments and exposure times, 
as in Fig. [14] The results are averaged over eleven reahsations. 

up to z ~ 0.4 in the hard band (for an integration time of 30 ks), 
while XMM-Newton reaches z ~ 0.9. 

As expected, only the most massive clusters can be iden- 
tified at high redshifts. In particular we note that the detec- 
tion threshold increases approximately by a factor of 5 between 
z ~ O.l and z ~ I - The mass of detected halos also increases in 
the hard and hardest bands because only the most massive halos 
emit such energetic photons. Finally, we recall that multi-band 
filtering is pointless for X-ray analyses because the noise in the 
different bands is uncorrelated, hence a multi-band filter would 
reduce to a singl e-band filter 



Compared to lFinoguenov et al.l (l2007h . who used a wavelet 
scale-wise reconstruction of the image (for details see 



Vikhlinin et al.l il998). we find a higher contamination in our 
sample, due to the fact that unlike iFinoguenov et al.l OOOTI) we 
do not use any optical information, but the amount of X-ray 
detections with S /N > 3 is similar. 

6. Correlation between X-ray and multi-band SZ 
detections 

We now describe the properties of common detections in the 
X-ray and SZ maps. We only consider X-ray detections in the 
soft band because their number was highest, and coiTelate them 
with multi-band SZ detections. A major problem is caused by 
the different angular scales of the detections in the different 
maps. Objects detected in X-rays are much smaller and more 
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numerous. Thus, many of them may overlap with the same SZ 
detection according to the exposure time. For XMM-Newton 
with an integration time of 100 ks, we found that a single SZ 
detection on average overlaps with at least ten X-ray detections. 
To solve this problem we only considered the brightest X-ray 
detection contained in the portion of the map associated to the 
SZ detection. 

Figure [16] shows the percentage of true detections in the 
X-ray observations which have a SZ counterpart. We show the 
cumulative percentage as a function of halo mass (top panels), 
and a differential distribution as a function of redshift of the de- 
tected halos (bottom panels). The left and right panels refer to 
light-cones with limiting redshifts z = 1 and z -2, respectively. 
Different line styles coiTespond to different instruments in the 
X-ray observations, while different colors refer to different in- 
tegration times. The curves are normalised to the total number 
of clusters present in both X-ray and SZ catalogues. The mass 
distribution of the common detections reflects the mass sensi- 
tivity of the SZ multi-band filter (see Fig. [T3] ). while the red- 
shift selection reflects the dependence of the X-ray sensitivity 
on redshift, which drops towards high redshift. 

In the lower panels of Fig. [161 we notice that the maxi- 
mum fraction of common detections occurs approximately at 
z = 0.3, due to the X-ray sensitivity. Several X-ray detections 
will be masked by a single SZ detection and thus not seen. This 
implies that combining SZ with X-ray maps will result in a halo 
number strongly biased by the angular extent of the SZ detec- 
tion. A given SZ halo will mask X-ray clusters projected into 
its close neighbourhood, independently of their redshift. This 
can be problematic for studies of substructures because it is 
impossible to identify substructures in the SZ map, and it will 
become difficult to study clusters at intermediate redshifts. This 
will be addressed in more detail in future work. 
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Fig. 16. Top panels: cumulative mass distribution for halos de- 
tected in both X-rays and SZ maps. Results are shown for light 
cones extending to z = 1 (left panels) and z - 2 (right panels). 
Bottom panels: redshift distributions for the same halos. The 
solid line refers to Chandra, the dotted line to XMM-A^evvfon. 
The black and red curves refer to the two different integration 
times (30 and 100 ks, respectively). 




Fig. 17. Upper panels: coiTelation between the tSZ and X-ray 
signals {y- and x-axes, respectively). The results are compared 
to the distribution of common detections, and are shown as a 
function of their S/N ratios (blue for < S/N < 1, red for 
1 < < 2, green for 2 < < 3, and black for S /N > 3). 
The left and right panels refer to the S/N values obtained us- 
ing the multi-band SZ and the (soft-band) X-ray filter, respec- 
tively. Lower panels: correlation between the values in the fil- 
tered maps (y-axis) and the original ones (x-axis) for the multi- 
band SZ filter (left panel) and the (soft-band) X-ray filter (right 
panel). Again detections are binned using their S/N ratios. The 
diagonal line would coiTespond to a perfect agreement between 
original and filtered maps. All plots refer to results for light- 
cones extending up to z = 1 . 

The upper panels of Fig. [17] show the correlation between 
the X-ray and tSZ signals of the common real detections 
grouped according to their S/N ratio significance, overlaid with 
the theoretical SZ and X-rays correlation as obtained in the 
original simulation. The analysis is limited to the soft band for 
XMM-Newton and an exposure time of 100 ks. 

As expected, detections with S/N > 3 occupy the 
upper right part of the contours (y > 10"^ and X > 
10"' ' erg/s/cm^/deg^) because they correspond to the most 
massive (and therefore brighter) halos. We also note that 
we do not find any detection with y < 10"^ and X < 
lO"'"' erg/s/cm^deg^ because these values of the flux corre- 
spond to objects whose mass is of the order of the minimum 
mass that the filters can detect. 

The lower panels show, the correlation between the filtered 
maps (y-axis) and the original maps (x-axis). On the left, we 
show the SZ multi-frequency filter, on the right the X-ray maps. 
We notice a substantial scatter for low Compton-y parameters 
in the SZ compared to the X-ray maps, corresponding to de- 
tections with very low S/N ratio. We also note that, while the 
filtered X-ray maps show unbiased values for S /N > 3, this is 
not strictly true for the multi-band SZ filter: the estimated am- 
plitudes are biased towards lower values. This may be due to 
the fact that the y6-model does not exactly reproduce the signal 
distribution of the SZ-halos. Finally, we found that the signal- 
to-noise ratios of the common SZ and X-rays detections are 
uncorrelated. 
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7. Conclusions 

We have tested the performance of linear filtering in optimising 
and combining galaxy-cluster detections in X-ray and SZ ob- 
servations. Based on a large A^-body hydrodynamical cosmo- 
logical simulation, we have studied the halo selection function, 
the contamination of halo catalogues and the mass sensitivity 
of the filtering technique. Our work extends earlier studies by 
taking gas physics properly into account and treating SZ and 
X-ray objects simultaneously with similar techniques. 

We constructed light cones up to redshifts 1 and 2 and sim- 
ulated typical observations of the X-ray emission (modelled af- 
ter the XMM-Newton and Chandra observatories) and the SZ 
signal (adapted to ACT) produced by the matter contained in 
the light cones. Our set of 11 different realisations covers total 
areas of ~ 264 deg^ and ~ 106 deg^ for limiting redshifts z - I 
and z = 2, respectively. 

We constructed and used single- and multi-band matched 
filters for halo detections in the synthetic SZ maps, and a 
single-band filter for the X-ray detections. We identified halos 
as S/N peaks in the filtered images and compared them with the 
parent distribution of dark matter halos contained in the origi- 
nal cosmological simulation. 

We find that multi-band filtering considerably reduces the 
number of spurious SZ detections and increases the fraction of 
halos reliably identified in the light cones. For instance, 50% of 
the halos with masses are found when using a 

multi-band filter, compared to 30% with the single-band filter 
The mass sensitivity of the multi-band filter is approximately 
four times higher than for the single-band filter As expected, 
the sensitivity is virtually independent of redshift. 

The fraction of spurious detections in the X-ray maps is 
typically low, and almost constant (about 5%) for a large range 
of S/N ratios. On the whole, the X-ray catalogues are more 
complete than the SZ catalogues, even when multi-band SZ fil- 
ters are used. However, the redshift distribution of the halos de- 
tected through their X-ray emission is strongly peaked at low 
redshifts, and only massive halos can be detected at high red- 
shift. This emphasises the importance of using complementary 
methods. 

Weak -lensing halos dete cted with the optimal filter pro- 
posed by iMaturi et alJ (l2005h reach a similar mass sensitivity 
as SZ and (soft) X-ray detections at redshifts where the lens- 
ing efficiency is highest, but clusters can be detected up to high 
redshifts. The contamination of weak-lensing halo catalogues 
is comparable to the X-ray catalogues, while the SZ contami- 
nation is much lower Weak-lensing halo catalogues are simi- 
larl y complete as SZ catalogues. For more quantitative details, 
see 
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Catalogues of common SZ and X-ray detections show that 
their mass function closely follows that of the SZ detections 
(due to the mass sensitivity of the filter), while their redshift 
distribution is very similar to that of the X-ray detections be- 
cause X-ray clusters can be observed more reliably at low and 
intermediate redshifts. 

Finally, we have shown that the detections represent the 
simulated halos in a virtually unbiased manner (see lower pan- 
els in Fig.fTTll. 
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